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The Priedel oscillations in the presence of transport currents 
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We investigate the Friedel oscillations in a nanowire coupled to two macroscopic electrodes of 
different potentials. We show that the wave-length of the density oscillations monotonically increases 
with the bias voltage, whereas the amplitude and the spatial decay exponent of the oscillations 
remain intact. Using the nonequilibrium Keldysh Green functions, we derive an explicit formula 
that describes voltage dependence of the wave-length of the Friedel oscillations. 



I. INTRODUCTION 

Transport properties of nanosystems, e.g., nanowires 
or single molecules, have recently been receiving signifi- 
cant attention mainly due to their possible application in 
future electronic devices i ' 2 i 3 i 4 These properties strongly 
differ from those of macroscopic conductors. The most 
important obstacle in theoretical investigations of the 
transport phenomena originates from the coupling be- 
tween nanosystem and macroscopic leads. Because of 
this coupling, analysis of the electron correlations is more 
difficult than in the equilibrium case. 

In nanosystems the charge carriers are usually dis- 
tributed inhomogeneously. There exist several reasons 
for such an inhomogeneity: 

(i) First, it may originate from a spatial confinement £ 
In analogy to the case of a quantum well, one may ex- 
pect that due to a small size of nanosystems, electrons are 
inhomogeneously distributed. In particular, recent scan- 
ning tunneling spectroscopy has shown the presence of 
the electronic standing waves at the end of a single-wall 
carbon nanotube^ 

(ii) Additionally, in the transport phenomena it may 
originate from the applied voltage.^^ In this case the 
system properties are determined by the chemical poten- 
tials of the left and right electrodes. Different values of 
these potentials may lead to an inhomogeneous distribu- 
tion as well. 

(Hi) Similarly to the macroscopic case, inhomogeneous 
charge distribution in nanosystems should occur in the 
presence of impurities 3 

(iv) Nanowires or molecular wires represent quasi-one- 
dimensional conductors. Therefore, phenomena typical 
for low dimensional systems, e.g., charge density waves 
may occur as welL 10 ' 11 ! 12 ' 13 ! 14 Recently it has been shown 
that the charge density waves are strongly modified by 
the bias voltage^ Apart from the low-voltage regime, 
they are incommensurate and the corresponding wave 
vector decreases discontinuously with the increase of the 
bias voltage. 

In this paper we focus on the impurity-induced inho- 
mogeneities. It is known that an impurity in the electron 
gas produces local changes of the carrier concentration, 
known as the Friedel oscillations^ that asymptotically 
decay with the distance from the impurity. The most 



of recent theoretical investigations of the Friedel oscil- 
lations concerned the influence of the electronic correla- 
tions, that is of crucial importance in one-dimensional 
systems! 16 i 17 ' 18 i 19 ' 20 It has been shown that correlations 
suppress the decay of the density oscillations j 17 ' 19 It is 
interesting that these oscillations give information about 
the impurities^ as well as the electron-electron interac- 
tion in Luttinger liquid systems i 18 ' 20 In macroscopic sys- 
tems, the Friedel oscillations are closely related to the sin- 
gularity in the response function for wave-vectors close to 
2kp, where kp is the Fermi wave-vector. If the nanosys- 
tem is isolated (or more generally, is in equilibrium), kp 
is a well defined quantity. However, in the transport ex- 
periments the nanosystem is coupled to two macroscopic 
leads with different Fermi levels and the difference be- 
tween these Fermi energies increases with the bias volt- 
age. Therefore, the meaning of kp is ambiguous. Since 
the properties of a nanosystem are determined by the 
chemical potentials of the left and right electrodes, one 
may expect that the Friedel oscillations should depend on 
the voltage as well. In this paper we analyze this depen- 
dence using the formalism of the nonequilibrium Keldysh 
Green functions. In particular, we derive an explicit for- 
mula for the voltage dependence of the wave-length of 
the Friedel oscillations. 

The paper is organized as follows: In Section II we dis- 
cuss a microscopic model and details of calculations. Nu- 
merical results are presented in Section III. Approximate 
analytical formulas are derived in Section IV. The last 
section contains a discussion and concluding remarks. 



II. MODEL AND THE CALCULATIONS 
SCHEME 



We investigate a one-dimensional nanowire with its 
ends coupled to macroscopic leads. The system under 
consideration is described by the Hamiltonian 



H — H c \ + H n 



3— el; 



(i) 



where H e \, ff nano and £f nano _ c i describe leads, nanowire, 
and the coupling between the wire and leads, respectively. 
We assume that electrodes are described by the free elec- 



tron gas, with a wide energy band: 



k.tr.a 



(2) 



where [i a is the chemical potential and a G {L,R} indi- 
cates the left or right electrode, hl — Hr = eV, with 
V being the bias voltage. c^ aa creates an electron with 
momentum k and spin a in the electrode a. The Hamil- 
tonian of the nanosystem is given by 



-ffnano = — t }^ d\ a d jcr + Uj^Tllc 
(ij)cr a 



(3) 



Here, d\ a creates an electron with spin a at site i of the 
nanosystem, rii a = d\ a di a and U is the impurity poten- 
tial. We have assumed a single impurity localized at site 
I. The coupling between the nanowire and the leads is 
given by: 



nano— cl 



k,i,a,<7 



(ffk^acj^C^ + H.' 



(4) 



In the following, we assume that the matrix elements 
5k, i, a are nonzero only for the edge atoms of the 
nanowire. 

The electron distribution has been determined with 
the help of the nonequilibrium Keldysh Green functions. 
Here, we follow the procedure used by Kostyrko and 
Bulka in Ref. 0- In particular, the local carrier density 
is expressed by the lesser Green function, 



i 



(5) 



This quantity, in turn, is determined by the retarded and 
advanced Green functions in the following way: 



E 

ae{L,R} 



G r (u)t a (u>)G a (w)f a (uj), (6) 



where 



f a (w) 



27r E 5 k,i,aSk,j,a%-£k,a), ( 7 ) 



and f a (uj) stands for the Fermi distribution function of 
the electrode a. The retarded Green function can be 
calculated from the following formula: 

d r (u)=\ul-H-t r (u)]~ 1 , (8) 



where H consists of the matrix elements of H aano , 
i.e., H 7nn t3m.n±l + U5 m ,i5 n: i and the retarded 

self-energy is determined by the coupling between the 
nanowire and the leads 



£» = 5 E 

ae{L,R} L 



Ip / dn^^- - if a (w) 

7T / U! — l2 



(9) 



0.3 



^ * ! ■ ' i! 



— eV = 





— eV = 


t 


-~eV = 


2t 




ii i 



40 



50 



60 70 
site 



80 



FIG. 1: (Color online) Occupation of sites in the vicinity of 
the impurity for a 129-site nanowire. The bias voltage is 
explicitly indicated in the legend. For the sake of clarity, the 
curves for eV = t and eV = 2t have been shifted downward 
by 0.1 and 0.2, respectively. 



III. NUMERICAL RESULTS 

We have solved numerically the system of Eqs. (5-9) 
for nanowires consisting of up to N = 129 lattice sites, 
with a single impurity in the middle of the wire. The only 
non- vanishing elements of T's have been assumed to be 



frequency independent 



ii 



r R (w) 



NN 



To, 



where the sites in the chain are enumerated from 1 to N. 
We have taken the nearest neighbor hopping integral t 
as an energy unit and assumed the coupling between the 
nanosystem and the leads as To = 0.1. The temperature 
of both the leads is fc^T = 0.01. Figure 1 shows the spa- 
tial distribution of electrons in the nanowire for U = 2 
and various values of the bias voltage. One can see strong 
oscillations in the vicinity of the impurity. However, due 
to the coupling to the leads the electron distribution vis- 
ibly differs from the standard Friedel oscillations: 



n(x) = n + A cos (Qx + rj) /x s 



(10) 



where the wave- vector Q = 2kp and the parameters ij 
and S depend on the interaction. In our case, the ac- 
tual value of Q has been obtained from the fast Fourier 
transform of the electron distribution n(x). 

Figure 2 shows the voltage dependence of the vector 
Q. In the equilibrium case (V = 0) Q = tt (with the 
lattice constant a = 1), so the charge oscillations are 
commensurate with the lattice. Since in the half-hlled 
case hp = n/2, this wave-vector remains in agreement 
with the standard relation Q — 2kp. However, when 
the bias voltage is switched on, the situation changes 
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FIG. 2: (Color online) The wave-vector Q of the charge os- 
cillations calculated by means of the fast Fourier transform of 
the numerical solution of Eqs. (5-9) (solid line). The dashed 
line shows the fit given by Eq. 1111 The numerical results have 
been obtained for the same parameters as in Fig. 1 The dis- 
creteness of Q(V) comes from the finite number of the lattice 
sites. 



and the oscillations are in general no longer commensu- 
rate. Moreover, one can see a strong dependence of Q 
on the applied voltage. Q is a monotonically decreasing 
function of voltage and vanishes for a sufficiently large 
V. Similar situation occurs in the transport phenomena 
through one-dimensional charge density wave systems. 15 
Our numerical results indicate that the obtained Q(V) 
dependence can be very accurately described by a for- 
mula 



eV = 4icos(Q/2). 



(11) 



The surprising simplicity of Eq. (|lip is very suggestive. 
In the following Section we present an approximate an- 
alytical approach that explains such a form of Q(V). It 
is applicable for arbitrary tight-binding Hamiltonian of 
nonintcracting electrons and holds true in a wide range 
of the coupling strength r . The numerical results pre- 
sented in Figure 2 allow us to test the applied approxi- 
mations. 



IV. ANALYTICAL DISCUSSION 

In the equilibrium case, the eigenstates of an isolated 
systems with periodic boundary conditions (pbc) are built 
out of plane waves. The Friedel oscillations are related 
to the maximum in the response function defined as a 
retarded Green function: 



x(Q,uj) = -((p(Q)\p^Q))), 



(12) 



calculated for U = 0. Here, 

P(Q) = ^ exp(iQRi)dl a d ia 

= XX+q<A°-' ( 13 ) 



where the summation is carried out over all momenta k. 
In the following we demonstrate that this quantity helps 
one to explain the dependence Q — Q(V) also in the 
noncquilibrium case. 

When the nanosystem is connected to macroscopic 
leads, the pbc become inappropriate since they do not 
reflect the geometry of the experimental setup. Then, 
the choice of open boundary conditions (obc) seems to 
be more appropriate. For U = 0, the Hamiltonian © 
with obc can be diagonalized with the help of the uni- 
tary transformation^ 



dl = 



N ■ 



jj2 sin ( kR i) d L> 



(14) 



where the wave-vectors k take on the following values 

7T 27T NTT 



k = 



N + V N + l' 



N + l 



(15) 



The specific form of this transformation accounts for van- 
ishing of the one-electron wave functions at the edges of 
the nanosystem. In this representation one gets 



ka 



where 



-It cos(fc). 



(16) 



(17) 



Although, the dispersion relation is exactly the same as 
for pbc, the values of k belong to (0, 7t) instead of the 1st 
Brillouin zone (— it, 7r). One can apply the above trans- 
formation also to the remaining terms in the Hamiltonian 
(P) and repeat calculations presented in Sec. II. The re- 
sulting equations have the same structure as Eqs. (5-9) 
with the real space variables i replaced by the wave- 
vectors k. In the new representation, the Hamiltonian 
matrix H is diagonal, however, the matrices T a take on 
much more complicated form: 



r a (w) 



fc P N + 



- V ff a (w)l sin{kRi) sm(pRj). 

1 L J i j 



(18) 

In order to analyze the Friedel oscillations we investi- 
gate the correlation function given by Eq. (TP2")) with 



P(Q) = ^cos(QRi)d\ a d ia 

icr 

= ^2 B Q( k >p) d L d P<?> 



(19) 
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where 

&Q{k,p) = - (Sp,k~Q - $p,Q-k + $p,k+Q ~ 8 p 27r -(k+Q)) ■ 

(20) 

Equations of motion allow one to calculate the correlation 
function that in the static limit takes on the form 

X (Q,w-0) = ]T ^fi(B* Q (p, q )(dld qa ) 

k,p,q,cr 

-B* Q (q,k)(dld pa ))+ X '. (21) 

The second term in the above equation, is propor- 
tional to g ((d'' c\d^ d)) and will be neglected in the follow- 
ing analysis. Such simplification is justified only for a 
weak coupling between the nanosystem and the leads. In 
order to demonstrate the validity of this approximation 
we have calculated numerically the resulting correlation 
function for nanowires consisting of 20 and 40 sites (see 
Fig. 3). The discreteness of the system is clearly visible 
for short nanowires, whereas for large systems the corre- 
lation function becomes smoother. In the latter case one 
can see that %(Q,w - ► 0) reaches its maximum value for 
Q given by Eq. (JTTJ , what justifies the applied approxi- 
mation x' — 0. 

The correlation function given by Eq. (|21[) is still too 
complicated for analytical discussion. The averages in 
Eq. (|2T|) are determined by the lesser Keldysh Green 
functions that, in turn, depend on the retarded ones. 
Because of the coupling between the nanosystem and the 
leads, off-diagonal elements of the retarded Green func- 
tion are nonzero [see Eqs. ((9]) and (|18j) ]. In order to 
proceed with the discussion of the Friedel oscillations we 
apply an additional approximation: 

(f a ) ~ 5 kp V (22) 

V / kp 

Then, all matrices in Eq. © become diagonal. In the 
next step we assume V —> 0, what allows one to calculate 
the integral over frequencies in Eq. ([5]). The resulting 
correlation function can be expressed as a sum of two 
Lindhard functions: 




FIG. 3: (Color online) Correlation function \ — x(Qi V) (Eq. 
I21[) determined numerically for 20-site (upper panel) and 40- 
site (lower panel) chain. 



x(Q^^Q) = x L (Q) + x R (.Q), (23) 

where 

X L W(Q) = £ \B Q (k,p)f fL(R) ^-_ fHR) (£p) . 

The Fermi distribution functions of the left and right 
electrodes read 

/l (efe) = / (ek~ , Ir (efe) = / {^k + , 

with f(x) = [exp(x/ksT) + 1] . The maximum of the 
response function occurs for such Q, that both the Fermi 



functions in the numerator in Eq. (|24p vanish simultane- 
ously. It is easy to check that for both x L (Q) an d X R (Q) 
this requirement is equivalent to Eq. pip . 

At this stage a comment on the approximation given 
by Eq. (|22[) is necessary. It is a crude and generally in- 
appropriate approximation that strongly affects most of 
the system's properties. In particular, it would strongly 
modify the current-voltage characteristics. However, the 
correlation functions calculated from Eqs. (|2"Tj) and (|2"3"]) 
are almost indistinguishable, what a posteriori justifies 
the use of this approximation for the discussion of charge 
inhomogeneities. This surprising result gives some in- 
sight into the physical mechanism of the charge distribu- 
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tion in nanosystems in the presence of transport currents. 
This distribution seems to be independent of the details 
of the coupling between the nanosystem and the leads 
however, it is determined by the fact that nanosystem 
is connected to two macroscopic particle reservoirs with 
different chemical potentials. 

V. CONCLUDING REMARKS 

Using the noncquilibrium Keldysh Green functions we 
have investigated the Friedel oscillations in a nanowire 
coupled to two macroscopic electrodes. We have derived 
a simple formula for the correlation function that de- 
termines the wave vector Q of the oscillations. The ap- 
proximate analytical expression fits the numerical results 
obtained from the Fourier transform of the electron dis- 
tribution very accurately. Our analysis concerns nanosys- 
tems described by the tight-binding Hamiltonian with the 
nearest neighbor hopping. However, it can be straightfor- 



wardly extended to account for arbitrary hopping matrix 
elements. 

The above discussion of the Friedel oscillations focuses 
on the voltage dependence of the wave-vector Q. We 
have found that the envelope of the charge density os- 
cillations is almost bias-voltage independent. It means 
that the remaining parameters characterizing the Friedel 
oscillations, i.e., the amplitude A and the spatial decay 
exponent 5, are determined predominantly by the inter- 
nal properties of the nanowire, whereas the wave-length 
of the oscillations depends on the bias-voltage. We be- 
lieve that investigations of the Friedel oscillations in the 
transport phenomena should allow one to get insight into 
many important parameters of the experimental setup. 
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